High-resolution temporal constraints on the dynamics of dark energy 
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We use the recent type la supernova, cosmic microwave background and large-scale structure data 
to shed light on the temporal evolution of the dark energy equation of state w(z) out to redshift 
one. We constrain the most flexible parametrization of dark energy to date, and include the dark 
energy perturbations consistently throughout. Interpreting our results via the principal component 
analysis, we find no significant evidence for dynamical dark energy: the cosmological constant model 
is consistent with data everywhere between redshift zero and one at 95% C.L. 



Introduction. It has been recognized for some time 
now that accurate reconstruction of the expansion his- 
tory of the universe is a crucial step towards understand- 
ing the physical mechanism behind the accelerating uni- 
verse. Early work on constraining the expansion history 
had concentrated on parametrizing the equation of state 
(ratio of pressure to density) of dark energy w(z) via one 
or two parameters and measuring them together with 
the energy density relative to critical SIde [1~5]. More 
recently, the program of reconstructing the expansion his- 
tory had been generalized to adding more parameters [6] 
and decorrelating them [7-9]. Nevertheless, the data at 
the time only allowed only up to three or four band- 
powers of w{z) to be considered, leading to "jagged" ex- 
pansion history, especially in the range 0.5 < z < 1. In 
the future, a number of "principal components" of dark 
energy — parameters forming a natural basis in which 
the function w{z) can be expanded — will be measured 
[10, 11]. 

In this paper, we utilize a new variant of the nearly 
model-independent approach to reconstruct w(z) of dark 
energy with the latest astronomical observations includ- 
ing type la supernovae (SNe la), power spectra of the 
cosmic microwave background (CMB) anisotropics and 
galaxy distribution from the large-scale structure (LSS). 
We perform a full likelihood analysis using the Markov 
chain Monte Carlo approach [12], and make sure to prop- 
erly take into account the dark energy perturbations [13]. 
Compared to previous work, our parametrization of the 
expansion history is more flexible and, as we argue, ro- 
bust and easy to implement. 

Method and Data. We consider the following cosmo- 
logical parameter set 



{w l , u b , oj c , 6 S , r, n s , log[10 10 A s ]} 



(1) 



where u>b = flbh 2 and uj c = Q, c h 2 are the physical baryon 
and cold dark matter densities relative to critical, O s is 
100 x the ratio of the sound horizon to the angular di- 
ameter distance at decoupling, r is the optical depth to 
re-ionization, A s and n s are the amplitude of the primor- 
dial spectrum and the spectral index respectively. 

The remaining parameters, Wi(i — 1,2, ..n), are the 
equation of state values at n specific fitting nodes {zi} 
— redshifts at which we vary these parameters indepen- 



dently. We use a cubic spline to interpolate between 
these nodes and obtain the full w(z) in the redshift range 
[0, z max ]. We set z max = 1.0 since the number and qual- 
ity of SN la data, and hence constraints on dark energy, 
dramatically weaken beyond this redshift. At z > z max , 
we assume w(z) = —1. We have also tried allowing 
w(z max < z < zcmb) to be a single new parameter 
free to vary; we found that the Markov chains do not 
converge because w(z > z max ) is difficult to constrain 
using current data. Comparing this approach to setting 
this parameter as — 1 , we found that the best % 2 remains 
nearly unchanged. Thus setting w{z > z max ) = —1 is 
safe in the calculation. Overall, our parametrization of 
w(z) can be summarized as: 



w{z) 




z > z n 

Z < Z u 
Z < Z-n 



,z G {zi}; 
,z £ {Zi}. 



(2) 



For the basic results, we choose n = 6 equation of state 
parameters (we discuss later the sensitivity to varying n). 

It is important to take account of dark energy per- 
turbations in the analysis, otherwise the results may be 
biased [3, 14]. In the conformal Newtonian gauge, one 
can derive the following perturbation equations [15]: 



-(1 +w){6- 3$) - 3H(c 2 - w)5 



-n(l-3w)9 



1 



fc 2 ( 



1 



(3) 
*) (4) 



where the over-dot represents the derivative with respect 
to conformal time, c 2 = SP/5p is the sound speed, W 
is the Newtonian potential, TL is the (conformal time) 
Hubble parameter, 8 is the density perturbation and 
the velocity perturbation. One easily sees that the per- 
turbations 5 and are divergent when w{z) crosses — 1 
even if nothing in the physical model should diverge. One 
way to solve this problem is based on the quintom model 
[16] which we adopt here. In particular we introduce a 
small positive constant e to divide the full range of the 
allowed values of w into three regions: I) w > — 1 + e; II) 
— 1 + e > w > — 1 — e; and III) w < — 1 — e. Neglecting 
the entropy perturbation contributions, for the regions I) 
and III) the equation of state does not cross —1 and the 
perturbations are well defined by solving Eq. (3) and (4). 
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FIG. 1: Reconstruction of w(z) using the datasets (I) (top panel) and (II) (bottom panel); see text for details. The red and 
black error bars correspond to 68 and 95% error bars respectively. ACDM model is shown with the blue dashed line. In the 
bottom panel, we also plot the reconstruction result with an alternative placement of fitting nodes (red points and curve; see 
the text for details). 



For the region II), the density perturbation <5, the veloc- 
ity perturbation 8, and their derivatives are finite and 
continuous for the realistic quintom dark energy models; 
therefore we set matching conditions in region (II), S = 
and 9 = 0. This is an approximate method to calculate 
DE perturbation during crossing regime without intro- 
ducing more parameters 1 . The error in this approxima- 
tion is controllable and we have tested it in our numerical 
calculations; with e ~ 10~ 5 we find that our method is a 
very good approximation to the two-field quintom dark 
energy model. For more details of this method we refer 
the reader to Refs. [3, 13]. 

In our calculations we take the total likelihood as the 
products of the separate likelihoods (Li) of CMB, LSS 
and SN la, i.e. defining xf = — 2 log A, we use 

2 _ 2 , 2 , 2 fr\ 

Xtotal — XCMB + XLSS + XSNIa ■ l°J 



Rcf[17] has proved that in FRW cosmology, the equation of state 
of a single fluid or a single scalar field cannot cross —1. To realize 
such crossing, one has to introduce at least one more degree of 
freedom. 



For the CMB data, we use the three-year WMAP 
(WMAP3) temperature-temperature and temperature- 
polarization power spectrum with the routine for com- 
puting the likelihood supplied by the WMAP team [18]. 
The LSS information we use consists of the galaxy power 
spectrum from the Sloan Digital Sky Survey (SDSS-gal; 
[19]); the luminous red galaxy power spectrum, also from 
the SDSS (SDSS-lrg; [20]) and the galaxy power spec- 
trum from the 2dF Galaxy Redshift Survey (2dFGRS- 
gal; [21]). To be conservative, for all galaxy data we 
have used only the information at k < fJ.lft.Mpc~ 1 , cor- 
responding to the linear regime. For SN la data, we use 
two datasets: Riess sample of 182 SNe [9] and ESSENCE 
sample of 192 SNe [22, 23]. We consider two dataset com- 
binations: 

(I) Riess-182+WMAP3+SDSS-gal+2dFGRS-gal; 

(II) ESSENCE-192+WMAP3+SDSS-lrg+2dFGRS-gal. 

In addition, for both datasets we also use the fol- 
lowing information: Hubble Key Project measurement 
H = 72 ± 8 km s" 1 Mpc _1 (l<r) [24]; baryon density 
information from the Big Bang Nucleosynthesis h 2 = 



3 



0.022 ± 0.002 (la) [25]; and a top-hat prior on the age of 
the universe, 10 Gyr < t a < 20 Gyr. 

For the MCMC calculation, we run eight independent 
chains each originally with O(10 5 ) elements, then thinned 
by a factor of 4. The average acceptance rate is about 
50%. We ensure the convergence of the chains by Gelman 
and Rubin criteria [26] and find R— 1 <~ 0(1O~ 2 ) which is 
more conservative than the recommended value R — 1 ~ 

©(io- 1 ). 

Results. We first study how the results change as we 
vary the number of and location of fitting nodes describ- 
ing w(z). We start out by placing one fitting node each 
at z = and z — z max = 1, then increase n, adding fit- 
ting nodes so that they are uniformly spaced in redshift; 
Zi = (i — l)/(n — 1), (n > 2, n e N). For each n and 
using the dataset (II), we calculate the improvement in 
X 2 relative to the ACDM model, A^ 2 = x 2 lin ~ Xacdm- 
As expected, we find that n = 2 already has a better 
fit than the ACDM model, and the fit improves as n is 
increased. When n is greater than 6, however, A% 2 flat- 
tens off (see panel (F) of Fig. 2). We therefore conclude 
that the data we use have enough power to constrain 6 
equation of state parameters. Therefore, at least when 
n < 8 is considered, n = 6 is the optimal choice since a 
higher n docs not improve the fit 2 . 

We also test the prescription of how to place the fit- 
ting nodes in redshift in the lower panel of Fig. 1 . The 
black points with error bars are for placement Zi G 
{0,0.1,0.3,0.5,0.7,1.0}; the red points with red curve 
are for placement z, G {0,0.2,0.4,0.6,0.8,1.0}. We see 
that the red curve almost coincides with the black points, 
namely, the dependence of result on placement of fitting 
nodes is weak. The above test and analysis has been also 
done with dataset (I) and a similar result is found. In 
what follows, we adopt the former, slightly non-uniform 
placement in redshift. 

Given the above two cosmological datasets, we explore 
the 12-D parameter space from Eq. (1). After marginal- 
izing over other cosmological parameters, we obtain the 
constraint of DE parameters w; as plotted in Fig. 1. In 
the upper panel, we find that almost all the mean val- 
ues of points are consistent with the prediction of ACDM 
model except for the 1st and the 4th points, at redshifts 
and 0.5. The equation of state at the present epoch 
is slightly favored to be below —1, while the equation of 
state at redshift 0.5 is greater than —1 at 68% C.L., show- 
ing a small "bump" . This bump most likely stems from 
a feature in the Hubble diagram of the Riess-182 dataset 
(I). To see this, we replace the Riess-182 and SDSS-gal 
data by ESSENCE-192 and SDSS-lrg data respectively 
and get the result in lower panel of Fig. 1. With this 
new dataset, the bump at redshift 0.5 disappears and 



2 We are aware of that the fit must further improve for some higher 
n because any dataset can be perfectly fit by some suitable, but 
possibly highly oscillatory, w(z). 



w(0.1 < z < 0.9) is consistent with -1 at 68% C.L.. The 
equation of state at redshifts and 1 deviates from —1, 
but not significantly. 

The equation of state parameters Wi are correlated, 
however, somewhat complicating their interpretation; for 
example, the highest correlation coefficient is —0.9, be- 
tween u>4 and u>5. It is therefore useful to de-correlate 
the Wi, as done in Ref. [7]. 

For each MCMC run, we compute the covariance ma- 
trix of the equation-of-state parameters, C = (wi — 
(wi))(wj - (wj)) T = (pp T ), using CosmoMC [12]. We 
then diagonalize the Fisher matrix F = C _1 , so that 
F = O t D O; here O is the resulting orthogonal matrix 
and D is diagonal. 

One can now rotate the parameters into a basis where 
the new parameters q are uncorrelated; q = Wp. There 
are many ways to do so [27] . Here we make two alterna- 
tive choices. 

Choice 1: Wi = F 1/2 : then 

C qi = (q 1 qf)=W 1 (pp T )Wf = 1. 

Choice 2: W 2 = O; then 

C q2 = <q 2 qj) = W 2 (pp T ) W T 2 = D 1 . 

Here Wi corresponds to the localized principal compo- 
nent (heretofore LPC) decomposition of the parameters 
p — the weights (rows of WJ are almost positive definite 
and fairly well localized in redshift [7]. Wi is usually re- 
scaled so that its rows sum up to unity, X^j=i(Wi)ij = 1, 
and we follow this practice. Note, however, that the phys- 
ical inferences do not depend on the normalization of Wi 
(or W 2 ) as both the parameters q and their values for a 
particular theoretical model change consistently. For ex- 
ample, we have tried a different normalization of Wi and 
explicitly checked that the parameters q and their values 
corresponding to the w(z) = —1 model change so that 
the hypothesis test of whether the measured parameters 
are consistent with w(z) = — 1 returns the same result. 

We find that all LPC band powers are consistent with 
the cosmological constant value at 95% C.L. Two out of 
six points deviate from ACDM at greater than 68% C.L, 
which is not statistically unexpected. This conclusion is 
easy to read off in the panel (A) of Fig. 2 where the q's 
are uncorrelated, and is independent of the choice of the 
decorrelating weights. In panel (C), we replace q x with 
q 2 (with W 2 normalized as Wi above), and again find 
consistency with the cosmological constant scenario. 

Summary and Discussion. Using the latest astronom- 
ical data of SNe la, CMB and LSS, we used the MCMC 
machinery to reconstruct the expansion history of the 
universe, and the evolution of dark energy, in a nearly 
model-independent fashion. We used the cubic spline in- 
terpolation between n values of the equation state Wj, 
fitting the expansion history more flexibly than essen- 
tially any description used on existing data so far. We 
found that n = 6 values of the equation of state in the 
range < z < 1 can be usefully constrained, and that 
increasing their number slightly does not improve the fit. 
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FIG. 2: Constraints on uncorrelated band powers of the equation of state w(z), using the dataset II. Panels (A) and (B) show 
the median values and the 68% and 95% constraints on the localized principal components q x , and their weights in redshift 
Wi, respectively. Panels (C) and (D) show the same for the principal components q 2 and their weights W2 (see the text for 
definitions). Panel (E) shows the likelihood distributions for the parameters q 1; while panel (F) shows the reduction in % 2 
relative to ACDM as the number of fitting nodes (for parameters q x ) is increased. 



Some comments can be made about our approach. We 
think the cubic spline leads a largely unbiased reconstruc- 
tion of w(z). An exception to this would be a highly 
oscillatory fiducial w[z), which however cannot be ro- 
bustly constrained with other approaches (and current 
data) either. Moreover, our results are not in conflict 
with Ref. [28] who claimed that only 2-3 equation-of-state 
parameters can be measured even from future surveys to 
better than about 10%, since we do not measure the six 
Wi to nearly such a good accuracy; see Fig. 1. Finally, 



one could certainly apply the same approach to the re- 
construction of the dark energy density Pde(z) as in [29], 
although we would argue that probing dynamics of dark 
energy requires specifically w{z). 

Our results do not show any significant evidence for 
the evolution of the equation of state with time, and are 
fully consistent with the cosmological constant scenario 
in the interval < z < 1 at the 95% C.L. We would 
particularly like to be able to test the dynamics of dark 
energy; for example, the "freezing" and "thawing" scalar 
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field models [30] which have different physical behavior 
and a different sign of dw/dz; however the accuracy we 
have right now is not sufficient to distinguish between 
these models [31]. Our method is straightforward to im- 
plement, and will produce sharp tests of the dynamics of 
dark energy once we have data from the next-generation 
surveys, such as Dark Energy Survey [32], Joint Dark 
Energy Mission [33], Large Synoptic Sky Telescope [34], 
and Pan-STARRs [35]. 
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